Emergence of superconducting textures in two dimensions 
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Self-organized regular patterns are ubiquitous in nature, and one of their most celebrated man- 
ifestations is the Abrikosov vortex lattice-: under an applied magnetic field, the homogeneous 
superconductivity becomes unstable and cast itself into a regular texture of the "normal" filaments, 
called Abrikosov vortices, immersed into a superconducting matrix. Its prediction and the exper- 
imental discovery became a breakthrough in our understanding of superconductivity and founded 
a new direction in physics. Here we show that the interplay between the superconducting order 
parameter and elastic fielda™™™*, which are intimately connected to the very existence of the su- 
perconductivity itself^, can result in a novel superconducting state dual to the Abrikosov state: a 
regular texture of superconducting islands. The fact that both patterns emerge within the framework 
of the Ginzburg-Landau description of superconductivity^ indicates that the formation of regular 
structures may be a generic feature of any phase transition. Emergence of superconducting island 
arrays is not specific to the effect of the elastic forces, but can be caused by any inherent mechanism 
generating long-range non-local interactions in the Ginzburg-Landau functional, for example, by 
the Coulomb forces. In particular, our findings suggest the formation of a superconducting island 
textures as a scenario for a superconductor-to-insulator transition in thin films. 



Sixty years ago the volume change accompanying the 
transition of a superconductor from the normal to the 
superconducting state was first observed^. This dis- 
covery - followed by finding the dependence of the su- 
perconducting critical temperature, T c , on the isotopic 
masses and the change in elastic constants at the tran- 
sition into the superconducting state^£ - geared the 
line of research crowned eventually by the triumph of the 
Bardcen-Cooper-Schricffcr (BCS) microscopic theory of 
superconductivity, which demonstrated an intimate con- 
nection between superconductivity and elastic properties 
of the material^. According to the BCS theory, T c is 
related to elasticity via£ 

k B T c = 1.13fiu D cxp (^^p^j i (!) 

where uj u is the Debye frequency, v is the density of 
states at the Fermi level, and V is the effective inter- 
action strength between electrons mediated by electron- 
phonon coupling. All three parameters, oj-d^vq, and V 
in ([I]) depend on pressure and comprise several micro- 
scopic effects including both the changes in the electron 
and phonon spectrum, as well as structural transforma- 
tions. The link between elastic properties and supercon- 
ductivity has been illustrated by numerous experiments 
revealing the influence of an external pressure on T c (see 
refs 11 and 12 for extensive reviews). 

Mechanical stresses, which affect superconductivity, 
are inherent to thin films attached to rigid substrates^ 
(Fig. 1). Indeed, enforcing on the film its own lattice 
spacing, the substrate impels internal strains in the film. 
On top of that, the rigid coupling between the film and 
the substrate can give rise to a peculiar scenario of the 
superconducting transition itself. Upon cooling the film, 




FIG. 1: The system. A sketch of a "soft" superconducting 
film deposited on a rigid substrate. Mechanical stresses in- 
duced by the substrate on the film due to mismatch between 
the lattice constants of the film and of the substrate, give rise 
to inhomogeneous superconducting state which emerges in a 
form of a regular array of separated superconducting islands. 



the superconductivity nucleates first in the regions where 
T c is elevated by fluctuations. Accordingly, mechanical 
properties, and, in particular, the lattice spacing in these 
regions should have changed. The rigid substrate, how- 
ever, obstructs expansions (contractions) of the film as- 
sociated with the local onset of superconductivity. As a 
result, additional local stresses emerge, which, in their 
turn, promote further growth (and the appearance of 
new) superconducting nuclei. The role of the substrate 
is thus twofold: first, it exerts elastic forces in regions 
where superconducting droplets appear first and, second, 
the substrate mediates the elastic coupling between re- 
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mote parts of the film, effectively transforming local dis- 
tortions into long-range elastic coupling. Appearing and 
evolving initial superconducting nuclei induce additional 
elastic distortions. As a result of this positive feedback 
between the superconducting order parameter and me- 
chanical stress, a spatial instability of the superconduct- 
ing state develops forming eventually a periodic pattern 
of superconducting islands (see Fig. 1). 

Our starting point is the time-dependent Ginzburg- 
Landau (TDGL) equatio n 8 ' 14 ? 15 combined with elastic 
stress balance equations (see Methods). This model of- 
fers a comprehensive description of formation and tem- 
poral evolution of the superconducting nuclei towards 
the final equilibrium configuration of the order param- 
eter. We choose a random initial configuration for the 
order parameter which mimics the random distribution 
of the superconducting nuclei expected to form in a 
homogeneously disordered film upon cooling it down 
to T c . However, the final equilibrium island pattern 
does not depend on the particular choice of the ini- 
tial configuration. The advantage of this approach is 
that all the richness and complexity of the microscopic 
elasticity-superconductivity interrelations is accounted 
for by the phenomenological coupling constant Uq = 
3a L K(dT c /dp), where p is the pressure, K is the bulk 
elastic modulus, and a L = (l/L)(dL/dT) is the linear 
expansion parameter (with L being the linear dimen- 
sion of the film). These quantities can be inferred from 
experiments and contain all the information about the 
electronic degrees of freedom and the lattice excitation 
spectrum^. 

In the absence of elastic interactions due to connec- 
tion to a substrate, one has the standard equilibrium 
solution of the conventional TDGL equation: a spatially 
uniform order parameter ip — tJjq describing the homo- 
geneous superconducting state at T < T c , and ip = 
at T > T c . The non-local elastic interaction, U(r — r'), 
coupling the values of the superconducting order param- 
eter at different points r and r' of the film [see Eq. (3) of 
Methods] distorts i/'o and can give rise to an instability 
of the uniform solution resulting in the formation of a 
regular island texture for certain values of the clastic pa- 
rameters. To investigate this process we solve the coupled 
TDGL and elasticity equations by numerical integration 
using a quasi-spectral technique. Figure 2 displays a se- 
quence of snapshots for the temporal development of the 
spatial structure of the amplitude of the order parame- 
ter. The time is measured in the units of the Ginzburg- 
Landau time t gl , with t gl = Trh/[8k B (T c — T)]. Starting 
from the random order parameter configuration, the sys- 
tem evolves through three clearly distinguishable major 
stages: (i) initial amplification of small fluctuations and 
emergence of an amorphous structure of islands; (ii) ap- 
pearance of a polycrystalline configuration of well sepa- 
rated superconducting islands; and (iii) slow relaxation of 
polycrystalline structure to a regular island lattice. The 
evolution of the modulus of the order parameter is quan- 
tified by the correlation function C(t) = Af^ 1 J^k IV^I 2 ; 



where ipk denote the Fourier components of the order pa- 
rameter, and the normalization factor M is chosen such 
that C = 1 when the island texture is fully periodic. 

The first stage is relatively fast: establishing of an 
amorphous island pattern takes only about I0t gl , while 
achieving a polycrystalline structure requires a 10 times 
longer period. In this intermediate time scale C(t) evolves 
logarithmically towards the polycrystalline state. In the 
final stage, for t > t p , see Fig. 2, where the island poly- 
crystal relaxes towards the regular lattice, the correlation 
function of the modulus of the order parameter ceases to 
be an indicative characteristic quantity. At longest time- 
scales it is rather the temporal development of the spa- 
tial distribution of the phase of the order parameter that 
characterizes the evolution of the system. For the chosen 
material constants the final state is a superconducting 
state with the uniform phase distribution corresponding 
to a long-ranged phase coherence. Figure 3 shows that 
the macroscopic phase-coherent state (which appears as a 
uniformly colored frame) is achieved via the motion and 
recombination of vortex-antivortex pairs which are ini- 
tially present en mass in the system due to coalescence of 
the independent superconducting nuclei in the first stage 
of the order parameter evolution. The vortex-antivortex 
pairs are clearly displayed by the "phase-cut" lines that 
appear as sharp color jumps from black to white. In the 
same time interval one distinguishes a well pronounced 
hexagon substructure which shows that different islands 
have different phases. The global phase coherence and 
therefore the superconducting state get established by 
the time of order 10 3 r GL , which is by an order of magni- 
tude longer than the time for establishing a robust dis- 
tribution of the amplitude of the order parameter. 

It is important that, while being a general phe- 
nomenon, the formation of an island texture requires 
the elastic coupling to be strong enough. Namely, the 
texture appears under the condition that the coupling 
constant Uq exceeds some temperature dependent criti- 
cal value U C (T). The function U C (T) is determined by 
means of linear stability analysis of the TDGL. Writ- 
ing down the Fourier transform of the elastic interac- 
tion term U(r) as U(q) = UoK.(q), where K,(q) is the 
scale free elastic kernel in the momentum space and 7 
is the diffusion coefficient in the Ginzburg-Landau equa- 
tion (see Method), one finds that the island state forms 
if (U K M - 1)(1 — T) - 7Q&/9.38 > 0, where K M is 
the maximum value of lC(q) at q = q M . This gives 
U C (T) = K- 1 + (jql /9.38/C M )(l - T/T c )-\ The lower 
bound for the coupling constant at which the homoge- 
neous state becomes unstable and the island texture can 
appear is U c = (7^/9.38 + 1)//C M . The left panel of 
the Figure 4 shows the phase diagram in the reduced 
T/T c -Uq/U c coordinates for the chosen system parame- 
ters which are given in the Methods. The right panel of 
Fig. 4 shows the behavior of the amplitude of the order 
parameter |A| at the isotherm T = 0.8T C as function of 
the coupling constant Uq. The protocol for finding this 
dependence is as follows: We start with a coupling con- 
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FIG. 2: Formation of the island texture. The framing sequence of snapshots shows the temporal evolution of the spatial 
distribution of the amplitude of the superconducting order parameter. The time is measured in units of the Ginzburg-Landau 
time, which for temperature T = 0.8T C can be estimated as t G l — 10 -11 seconds for T c ~ IK. The first frame at time t = 
shows an initial random configuration of the order parameter. One can distinguish three stages of the evolution: (i) emergence 
of an amorphous structure of islands (for t < t a ~ 10); (ii) formation of a polycrystalline islands texture (for t a < t < t p ~ 100) 
and (iii) relaxation to a long-range ordered island lattice (for t > 10 4 ). The color bar for the amplitude of the order parameter is 
presented beneath the initial frame (note the change of the scale for t > 3). The upper part of the central panel shows the time 
evolution of the normalized order parameter correlation function, C(t) — Af^ 1 ]>^ fc IVM 2 ! vs i on a logarithmic scale (tpk denote 
the Fourier components of the order parameter and the normalization factor TV is chosen such that C = 1 when the island 
texture is fully periodic). At intermediate times between t a and t p the correlation function shows transient logarithmic behavior 
(highlighted by the straight line) and exhibits for t > t p a slow convergence to unity. The lower part shows a perspective view 
of the height profile of the amplitude of the order parameter corresponding to a small region of the perfect lattice appearing at 
the final stage of solution of the TDGL equation. The simulations were done with the coupling constant of Uo = 2.23[/ c and 
the thickness of the film of 0.8£, where £ is the superconducting coherence length at zero temperature. Defining an "island" 
as an area within which the amplitude of the order parameter exceeds half of its maximal value, we find their size to be about 
2.5£ and the distance between the centers of the islands to be 12£. 



stant Uq < U c and an initial random distribution of the 
order parameter, and let the system evolve until a sta- 
tionary distribution of A is achieved. After that Uq is in- 
creased and the system evolves again until a new station- 
ary state is achieved (adiabatic increase), and so forth. 



Upon crossing the phase boundary a bifurcation in the 
| A (Uo) | -dependence corresponding to the islands forma- 
tion occurs and a finite difference between the maximal 
and minimal values of the amplitude of the order param- 
eter appears. In other words |A(£/o)| shows finite jumps 
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FIG. 3: Temporal evolution of the amplitude and phase of the order parameter. Four sequential snapshots of the 
amplitude and phase of the order parameter taken at times 100, 200, 300, and 1000 measured in units of tgl are presented (the 
same values of the system parameters as in Figure 2 were used). The first three frames correspond to nonequilibrium states 
characterized by an inhomogeneous phase distribution; in this stage the formation of a regular texture of the order parameter 
amplitude goes along with the recombination of vortex-antivortex pairs. The first frame displays four sets of vortex-antivortex 
pairs with the endpoints A-B, C-D, E-F, and G-H. In the corresponding phase frames these endpoints confine the line of 27T 
phase jumps visible as sharp color change from white to black (the periodic boundary conditions were used, such that the line 
G-H goes over the upper- and lower edges of the frame). The phase cut lines A-B and E-F disappear at t — 200, and at t = 300 
only the G-H line remains. In the nonequilibrium states, the hexagon-like structure of the phase distribution reflecting that 
on the way to equilibrium the phase at the islands are different is clearly distinguishable (each hexagon represents the phase 
of an island). At t — 1000 the phase becomes homogeneous across the sample and a global phase-coherent superconducting 
state establishes. However, for a different set of parameters the phase equilibration might be very slow or even stopped in the 
presence of impurities or dissipation, resulting in a random pattern of phases for each island, shown in the inset of the last 
phase- frame. 



to the |A m i n (t/ )| and |A max (E/o)| branches at some Uq 
appreciably beyond the phase boundary On the descend- 
ing part of the cycle, \A min (U )\ and |A max (t/ )| merge 
at smaller value of Uq closer to the phase boundary, i.e. 
at smaller Uq. The change of Uo is indicated by the ver- 
tical arrow in the left panel of Fig. 4. Shown in the right 
panel inset is the hysteretic behavior of the order parame- 
ter corresponding to sweeping the temperature across the 
phase transition line. The protocol is nearly the same as 
for the coupling constant variation scheme; the temper- 
ature is first decreased from T = T c down to T = 0.8T C 
and then increased back at the fixed coupling constant 
Uq = 2.6U C . The temperature sweeping is indicated by 
the horizontal arrow in the left panel of Fig. 4. Note 
that the island structure persists even as the system is 
heated above T c ; this is the manifestation of the positive 
feedback between the superconductivity and elasticity in 
this films bonded to the substrate giving rise to a local in- 
crease in superconducting transition temperature within 
the islands. The observed hysteretic behaviors indicate 
that the transition between the homogeneous supercon- 
ducting and the island texture states is of the first order. 

The pattern formation and, in particular, the spon- 
taneous emergence of electronic nanometer-scale struc- 
tures, due to the existence of competing states is ubiqui- 



tous in nature and is found in a wealth of complex sys- 
tems and physical phenomena ranging from magnetism, 
superconductivity and superfluidity to liquid crystals and 
biological dynamic o 15 ! 17 . The Ginzburg-Landau equa- 
tion, one of the most universal models in theoretical 
physics, offers an adequate quantitative description of 
the pattern formation. Here we have revealed the in- 
stability and emergence of an electronic textures, due to 
coupling of the superconducting order parameter with 
clastic stress fields, generating a non-local term in the 
Ginzburg-Landau functional. It is noteworthy that the 
general form of the nonlocality in the Ginzburg-Landau 
equation, induced by the interplay between the local non- 
linearity and the nonlocal couplings, is universal and 
can result from various physical mechanisms. In par- 
ticular, while the results presented here are stable with 
respect to weak disorder, we expect that the interplay 
of strong disorder with Coulomb interactions will give 
an essential contribution to the nonlocal coupling in the 
Ginzburg-Landau functional. Note, in this connection 
the interesting discussion of the effects of disorder lead- 
ing to a spatially inhomogeneous distribution of the order 
paramete r 18 i 19 i 20 i 21 i 22 . We thus conclude that a consis- 
tent consideration of the superconductivity in thin films 
has to include the effects of the interplay between strong 
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FIG. 4: Phase diagram and hysteretic onset of the 
island texture. The left panel shows the phases in the tem- 
perature T - coupling constant, Uo, plane. Above T c the film 
is in the normal state, whereas below T c the film can be either 
a homogeneous superconductor or a textured superconductor 
consisting of isolated islands. The phase boundary is deter- 
mined by the stability condition of the TGLE (see text). The 
right panel shows the amplitude of the order parameter A 
as function of the coupling constant Uo- The onset of the 
island texture (crossing the phase boundary) is marked by 
the bifurcation point where A max and A m i n start to diverge. 
The observed hysteretic behavior - the island structure forms, 
upon adiabatically increasing the coupling, at larger value of 
Uo than the reverse transition, upon decreasing Uo, to the spa- 
tially homogeneous order parameter state - indicates that the 
formation of the island texture occurs via a first order tran- 
sition. The path of change of Uo is indicated by the vertical 
arrow in the left panel. The inset in the right panel shows the 
hysteretic effect corresponding sweeping temperature at fixed 
Uo = 2.6U C , first, from T c down to 0.8T C , and then heating 
up (the process shown by the horizontal dotted line in the 
left panel). It is noteworthy that the island structure persists 
even upon heating the system above T c of the film material, 
due to the increased local T c within the islands. The dot- 
ted line indicates the usual \J\ — T/T c behavior of the order 
parameter without elastic interactions. 



disorder with Coulomb interactions and elasticity. 

The observed island texture is an array of granules of 
a "good" superconductor immersed into a sea of a "bad" 
one with much smaller amplitude of the order parameter 
and is equivalent to an array of small superconducting 
droplets connected via weak links, i.e., a Joscphson junc- 
tion array (J J A). In our simulations, where thermal and 
quantum fluctuations are absent, the Joscphson coupling, 
however weak, eventually establishes a global phase co- 
herence and thus global superconductivity. The fluc- 
tuations, dissipation, and disorder can break down this 
phase uniformity: the intermediate stages of phase dis- 
tributions of our simulations represent the final states of 
the "real world" systems. Thus the possible eventual con- 
figuration of the order parameter phase (see the "Gaudi- 
mosaic" phase distribution in the inset of the last phase- 



frame of the Fig. 3) can correspond even the insulating 
state. This is in accord with the general result that JJAs 
turn insulators if the Joscphson coupling between the su- 
perconducting islands becomes sufficiently weak^ 3 -. Thus 
our finding that the transition into a superconducting 
state in two dimensional superconductors elastically tied 
to a substrate can occur via the formation of nanoscale 
superconducting island pattern offers a possible scenario 
for the superconductor-to-insulator transition observed 
in various thin superconducting fihn a 24 i 25 i 26 i 27 i 28 i 29 . This 
supports the conjectur e 30 ' 31 that the formation of a reg- 
ular array of superconducting droplets is an inherent 
property of the critical region of the superconductor-to- 
insulator transition. 
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APPENDIX A: ADDENDUM 

1. Ginzburg-Landau equation. 

The two-dimensional time-dependent Ginzburg- 
Landau equation (TGLE) describing the behavior of the 
film is written, in dimcnsionless variables as: 

d t tp = ai; - f3\ip\ 2 ip + jV 2 4> - 5\i;\ 4 i; , (Al) 

where ip = ip(r,t) is the complex, dimcnsionless super- 
conducting order parameter normalized with respect to 
T c . We consider the static case where the external fields 
are absent, thus the electro-magnetic potentials can be 
absorbed into the order parameter. We include a small 
sixth order term in the GL free energy to ensure a numer- 
ical stability. The coefficient a = (r) is a functional of 
the order parameter and includes the effects of long-range 
clastic potentials, temperature T, and disorder: 

cty(r) =a [T c (r)-T]+ / dr'U(r - v')\i>(r', t)\ 2 , (A2) 

where T c (r) = T c (l + r(r)) describes quenched fluctua- 
tions (r(r) £ [—A; A]) of T c and U(r — r') is the non-local 
kernel generated by elasticity interactions. The tempera- 
ture T is measured in units of superconducting transition 
temperature T c , and the unit of the length is the super- 
conducting, zero temperature coherence length £. 
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2. Elastic interaction. 

To couple elasticity and superconductivity we utilize 
the pressure dependence of T c (p) which is taken in a lin- 
earized form as 



T c (p =Po + Ap) = T c (p ) + Ap 



d T c ( Po ) 
dp 



(A3) 



where Ap is the change in the internal pressure induced 
by elastic forces. The latter is derived from the condition 
of the continuity of the deformation fields u^- s ^ and the 
stress balance at the film-substrate interface (the z = 
plane) , 



<4 s) (0,0,0) 



aL p) (0,0,0) 



for 



x,y,z, 



and 



uf } (0,0,0) = uf> (0,0,0) for i = x,y,z, 

and at the film-vacuum (z = d, d is the thickness of the 
film) interface: 



o£>(o,o,d) = o 



for 



x,y,z. 



(A4) 



The superscripts (s) and (p) refer to the film and to the 
substrate respectively. 

These nine coupled equations are solved using a plain- 
wave Ansatz, yielding the exact expressions for the defor- 
mation fields which, in their turn, are used for calculating 
the internal pressure of the film. The Fourier transform 
of the pressure in the film then reads: 



P 



d y u y 



d z u ( z 



(A5) 



Using this expression together with the relation Ap = 
—KAV/V in Eq. (|A3|) one finds the parameter a of the 
TDGL in the form given in (|B7|) . We define the in- 
teraction potential U (r) through its Fourier transform 
U )C(q), where U = 3KAa L [dT c (p )/dp\ and K(q) is 
the scale free part of p^ (<l)/\ip q \ 2 [note, that p^(q) oc 
IV'gl 2 ]- Here K is the bulk modulus, the linear thermal 
expansivity of the film is Aa^, and the relative volume 
change due to the change in pressure Ap is AV/V . 



APPENDIX B: SUPPLEMENTARY 
CALCULATIONS 

1. Construction of the elastic interaction potential 

U 

Here we derive U (r) for the problem of a thin super- 
conducting elastic film coupled to a rigid substrate with 
different elastic properties and lattice mismatch. In the 
following, we use the superscripts (s) for all quantities 
related to the film, and (p) for the substrate. Both mate- 
rials are described by the three dimensional displacement 
fields /u^ . We start our consideration using a plain 
wave Ansatz for these fields: 



,(•) 



As) 



,(p) 



e w.*+»«»» (e q * z A hx + A 2<x e- q ** 



i {s) 
l 0,y 



-qzz) 



Q e qz z +iqx x +i-qyy 



where the misfit strain on the film is captured by the real 
exponential terms of the film in z-direction. Uq x and Uq y 
are special solutions of the two-dimensional problem due 
to superconductivity, which we will discuss later. In order 
to determine the nine free parameters Cj, A\ t i, and A2^ 
we need nine equations. First, we consider the interface 
film-vacuum at z = d, where d is the thickness of the 
film. The stress boundary condition for a flat surface 
demands: a xz = a yz = a zz = 0. So, for the film we get 
from Hooke's law 



T« = 

xz 

7< S > = 
>yz 

T W = 



M W (d z u x ^ + (l-2a )d x u^ 
V {s) (d z u^ + (l-2a )d y ui s) 



2/i 



1 + ^ s >) d z u^ + (d x u^ + 9j,4 S) ) 



3. Numerical simulation. 

The TGLE (|B9|) is solved on a fine discrete grid by the 
numerical integration in real- and Fourier spaces (quasi- 
spectral split-step method), taking into account the full 
non-local elastic potential. We use periodic boundary 
conditions and N = 512 2 grid points for a system of 
the size L 2 = (200£) 2 . The temperature is chosen as 
T = 0.8T C and the coupling constant of Uq/U c = 2.23 
(in Figs. 2 and 3) for a film of thickness d = 0.8£. 
The used GL parameters are ao = 4.37, (3 = 1/2, 
7 = 0.01, and 5 = 0.1. The parameters of the elas- 
tic kernel are p^ = 0.5, /^ p ' = 5.0 (shear moduli) and 
V ( S 'P) l(\ — 2v( s ' p ">) — 1.6 (modified Poisson numbers) - 
see supplementary materials for details. 



where the effect of a lattice mismatch of both materi- 
als is captured by the cro-term. In general, the lattice 
constants of both materials are different, i.e., the mis- 
fit parameter rj = (a^ — a^)/a^ is non-zero, where 
and a^ p ' are the lattice constants of the film and 
substrate, respectively. This lattice mismatch leads to 
a compression of the film-lattice at the interface, such 
that the lattice spacings match there. We note, that this 
description not only applies to crystalline structures but 
also for amorphous materials, where the "lattice" con- 
stants are averaged quantities. 

The coefficients fx^ s ' and v' 8 ' are the shear modulus 
and the modified Poisson number of the film [v^ = 
^^/(l — 2j>' s )), where is the Poisson number] . In the 
following we use the fact, that all interface equations are 
independent of the x and y coordinate and set x = y = 0. 
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Therefore, the first three of the nine equations are 

a$ (0,0,d) =0 for i=x,y,z. (Bl) 

Next, we consider the interface between film-substrate 
at z — 0, for which we also need the expressions for the 
stress tensor elements 



2^ 



rip) 



1 



(d x 



,(p) 



dyU^ 



At this interface we need to fulfill six more conditions, 
where three of them are resulting from the stress balance 

a\f (0, 0, 0) = a {p) (0, 0, 0) for i = x,y,z, 

and the last three equations are the continuity equations 
for the displacement fields 

u\ s} (0, 0,0) = u (p) (0, 0, 0) for i = x, y, z . 

These nine equations can now be solved in order obtain 
the free parameters of our Ansatz. 

The two special solutions for vS*> in x- and y-dircction, 

(s) 

u [ , already used in the Ansatz for the deformation fields 
of the film are obtained through solution of the extended 
stress-strain balance 



/3MAuW + v(VuW) = 7 (s) V|V(r)r , 



(B2) 



which takes into account that the systems energy is the 
sum of the superconducting energy and the mechanical 
deformation energy. This relation introduces two phe- 
nomenological parameters: 7^' and /?W which together 
define the coupling constant of the elastic interaction and 
superconductivity which is determined explicitly below. 
Eq. (|B2[) can be solved and gives 



1 0,i 



-lis 



f{l)li 



(1 + A)q 2 ' 



(B3) 



for i = x, y. 

Using the solutions for all coefficients we find the in- 
ternal pressure of the film 



p {s) (qx,q y ,qz) 



1 

3d 



dz 



(d x i 



dyUP 



d z ui s) 



(B4) 

This expression simplifies to p 1 • ^(q z ) if we use the homo- 
geneous parametrization q x = q z cos(8), q y = q z sm(6). 
Due to the special solution for u W the pressure p^ (q) is 
proportional to V'gl 2 : where ip q are the Fourier compo- 
nents of the order parameter. 

Now, we need to consider the connection of elasticity 
and superconductivity. For that, we first expand T c in 
changes of (internal) pressure: 



T c {p = po + Ap) = T C ( PQ ) + Ap 



dT c ( Po ) 
dp 



(B5) 



The pressure change Ap is related to a volume change 
AV as K = — Vdp/dV, where K is the bulk modulus. 
Therefore we can write 

Ap = -K^- = -3K—^- cx SKoil , (B6) 
V Lj 

where AL/L is the relative length change of the system 
and ctL the linear thermal expansivity of the film. 

Using the above result in the pressure-expanded criti- 
cal temperature, the parameter a in the TGLE obtains 
thus the non-local form shown in the Methods Sum- 
mary section: 

c^(r) =a [T c (r)- T} + J dv'U{r - v')\i>(r', t)\ 2 , (B7) 

where T c (r) includes also spatial variations of T c due to 
weak disorder. 

If we define the scale free clastic kernel 



K{q) 



(1 + 13^) p(«)( g ) 
7 ( s ) |V(r)I 2 



(B8) 



where the explicit form of the elastic kernel is quite 
involved and given in section IB 31 Finally, we write 
the Fourier transform of U(r) as UolC(q), where Uq = 
3KAa L [dT c (p )/dp} 



2. Numerical realization 

The TGLE equation 

d t ^ = aiP- P^^ + ^tp-S^Tp, (B9) 

is solved by a quasi-spectral split step method for periodic 
boundary conditions on a two-dimensional square grid 
with grid size of N 2 for a system size L 2 . Split step 
means, that for each time step the equation is solved 
partially in real space and partially in Fourier space (the 
diffusion part)^. In the first step we calculate 



4 h3 (t + At) = e^ a -^^\^{t)f 



(B10) 



where the non-local part of a is calculated by Fast- 
Fourier- Transformation (FFT) upfront. 

Then, in the second step, we first Fourier transform 
ipij(t + At), apply the diffusion kernel, and transform it 
back. This way we avoid mostly all complications of the 
diffusion equation^. 

In general the TGLE without long-range potential has 
only two stable solutions: a homogeneous one and a 
striped phase. Even if a long-range potential is present, 
the stability of these two solutions is only destroyed un- 
der certain conditions which can be found out by linear 
stability analysis (see next section). E.g., a Coulomb po- 
tential or even screened Coulomb potential cannot cre- 
ate this kind of instability. However, the elastic po- 
tential does, if the model parameters are chosen ap- 
propriately. Another possible non-local potential which 
destroys the homogeneous/striped solutions is a box- 
potential^ U(r) = U @(l - \r\/a). 
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TABLE I: List of parameters for the elastic model. 



a 




thickness of the film 


/' 


v) 


tiViPfir innrliiliiQ 

□ llLtll lllU^l LID 


<7(1 




deformation stress 


«/<* 


p) 


modified Poisson number 


Uo 




potential strength 


K 


p) 


coupling constant 




stress tensor 


4 


p) 


displacement fields 



Ai t i, A2,%, Ci 9 free variables for the elastic equations 

3. Elastic kernel and simulation parameters 

In order to determine whether the elastic long-range 
interactions lead to an instability of the homogeneous 
solution ipo of the GL equation, we need to check the 
linear stability of the equation, by expanding the order 
parameter around ipg. After this expansion and Fourier 
transformation, we find that the stability depends on the 
behavior of the function 

S(q) = (2U JC(q) - 1)(1 - T/T c ) - 79 2 /9.38 . (Bll) 

We know that IC(q = 0) = and JC(q — ► oo) = 1/3 and 
therefore S(q = 0) < (for T < T c ) and S(q) -> -oo for 



q — > oo. If S(q) changes its sign at intermediate q the 
homogeneous solution becomes unstable, see Fig. [5] 




FIG. 5: Linear stability. Plot of the stability finction S(q) 
with elastic parameters chosen as described in the text for 
different temperatures and Uo = 2.6U C (cf. inset of Fig. 4). 



For completeness we write the scale free elastic kernel 
K, explicitly, depending only on dimensionless material 
parameters: 



+d q \e Ad V^ (s) (2M (p) + 2m Cp) i' Cp) cto - At (a) i? (s) ) 

-2e 2d " (V?) 2 [l + 2*»] [l + 2i/«<7o] - V*> a - M W M W#W( V W + „M _ 2i»<r ))] 
+M (S) (e d "-if [V'V') (e 2d " - l) 0« [<7 - 1] 
+ fiCri [i + 2i/( s ) (2 - 2tf( s ' - 3cr ) + 2l» [2tf( a > + l] [a - 1] - 2e d 1 (i» hu&> [a - 1] + 2cr - 3) - l) 
- e 2d q (3i3 (s) +2 + 2i» [a - 1] (V s) [4<t - 2] - l))J] 
+4/i (p) V a >e 2d 9 [l + 2^( p) ] [4a - 3] sinh (d g)} / 

{3d f/ {V s)2 ( e 2d9 -l) 2 tf (s)2 

+2/>V s) (e 2d 9 - l) <? (s) (l + ^ (p> + v (a) ~ v {p) vo + e 2d q [l + i>V l) 

+M (p)2 (l + 2!») (e 4d +2l /( s ) [ CTo - 1] - 1 -2e M « [l + 2*»<r ] ) } } (B12) 

I 



with $( s ) = 2v^ [do — 1] — 1 and parameters described 
in table EU 

For the simulations (cf. evolution of amplitude and 
phase in Figs. 2 and 3) we used the typical values: 
T = 0.8T C , U = 2.2U C , 7 = 0.01, f3 = 0.5, d = 0.8£o, 



/i' s ) — 0.5, = 5 (the substrate is more rigid), 

;,(>) _ u (p) — 2.6, <To = (the lattice mismatch in the 
normal state is not relevant), L = 200£o (linear dimen- 
sion of the system), and N = 512 (number of discrete 
grid points per dimension). 
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